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positions and depths of these minima are regulated by the counterion density in 
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the coexisting phases for DNA bundles are calculated. A DNA-condensation and 
redissolution transition and a stable mesocrystal with an intermediate lattice constant 
for high counterion concentration are obtained. 
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1. Introduction 

Long DNA is naturally found in a dense form in most biological systems due to 
the presence of compacting agents. In vitro the most used agents are polyamines 
such as trivalent spermidine (Spd) and tetravalent spermine (Spe) p^. These agents 
play a key role in maintaining cellular DNA in a compact state [3 0], modulate ion 
channel activities of cells, are essential for normal cell growth [3] and can effectively 
be applied in gene delivery and in the field of genetic therapy. Under physiologic 
ionic and pH conditions, the polyamines are positively charged and hence DNA is 
their prime target of interaction. The molecular mechanism of polyamine function 
in DNA condensation is presumed to involve neutralization of the negatively charged 
DNA backbone by the positively charged amino groups of Spd and Spe |S] . Experimental 
results and counterion condensation theories indicate that non-specific interactions are 
predominantly electrostatic between polyamines and DNA phosphates [HI El IHl E] • Thus, 
the electrostatic shielding of phosphates facilitates close helix-helix surface contacts and, 
ultimately, DNA condensation through the correlation attraction jTUj, the attraction 
between strongly correlated counterion layers on the adjacent DNA surfaces. 

In the last decade, different experiments |31 CH 1121 UHl HH UHl HE] have shown 
evidence of DNA redissolution (i.e. DNA unbinding) with increasing concentration of 
polyamines and several of their structural analogues. First, addition of a certain amount 
of multivalent salt causes the collapse of single DNA or the bundling in multi-columnar 
structures. Upon addition of more salt, the polyelectrolyte redissolves and DNA assumes 
again an unbundled conformation. There are three experimentally well established 
features of the redissolution phenomena: i) a linear relationship is found between the 
threshold concentration of multivalent ions Cc, when the onset of DNA condensation 
takes place, and the initial DNA concentration co^a (the DNA concentration in solution 
free of multivalent ions); ii) the decondensation threshold Cd of multivalent ions, when 
the condensed DNA returns back to solution, is almost independent of monovalent salt 
concentration c^; iii) between the threshold values Cc and Cd the DNA fragments show 
two coexisting liquid crystalline phases. Concerning the first item, the onset of DNA 
condensation is usually explained by the correlation attraction between almost neutral 
structures. Thus, the precipitation induced by trivalent or tetravalent ions is not a 
consequence of the intrinsic structure and flexibility of DNA, but is a common feature 
of a polyelectrolyte solution. The threshold value Cc grows with increasing monovalent 
salt concentration j2l El El IHj • A mono- molecular DNA collapse into a neutral 
toroidal structure occurs in highly dilute solutions of long DNA chains jHl IH] , whereas a 
multi-molecular aggregation is generally observed in more concentrated DNA solutions, 
regardless of the DNA length [201 1^- The second item, the DNA redissolution at 
Cd, is currently under intensive debate in the colloidal community with different, and 
sometimes conflicting, explanations [121 CSIIIHIIIZI I22| • For instance, in Ref. [12] it is 
argued that after precipitation the increased binding of polyamines will make the DNA 
hydrophilic enough to dissolve in water. Other experiments [22] show that the DNA is 
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still in a condensed state when polyions are added beyond the threshold concentration 
Crf, but it is more finely dispersed in the solution. In Ref . ^Tj the reentrance is explained 
by resorting to DNA overcharging phenomena, which take place roughly in the middle 
of the condensation plateau. DNA is claimed to experience negative electrophoresis and 
move opposite to the conventional direction near the reentrance transition. However, 
experiments of Raspaud et al P do not confirm this claim. In Ref. [12] it was suggested 
that the redissolution is very sensitive to the relationship between the condensation of 
multivalent counterions on the polyelectrolyte and multivalent counterion-monovalent 
colon association (Bjerrum pairs). Thus, if the chemical potential of the multivalent 
counterions is low, they more likely adsorb on the DNA and overcharge it. On the other 
hand, if the chemical potential is high, the counterion-coion associations are more likely 
to appear in solution. A resulting condensation of Bjerrum pairs creates consecutive 
layers of oppositely charged ions around the DNA molecule [HI (231 • 

While the first two above mentioned items have been studied in considerable detail, 
much less attention has been paid to the third item, namely the coexistence of two 
different liquid-like structures in DNA condensates. In a series of experiments, Livolant 
and colleagues demonstrated that the spermidine and spermine ions are capable of 
provoking several liquid crystalline forms of fragmented DNA JT]. Similar results were 
published recently by Saminathan et al in Ref. [Sj- 

In this paper we investigate the condensation and redissolution of DNA on 
a molecular level by using computer simulations of the primitive-model electrolyte 
with explicit tetravalent counterions and monovalent salt ions. We trace back the 
condensation and redissolution to the distance-dependent effective potential U{R) 
between two parallel DNA molecules, where R is the radial distance between their 
two centers. Using two-dimensional liquid-state theory for the fluid and lattice sums 
for the solid phases, we calculate the phase diagram for columnar DNA assemblies. A 
previous account of the results was already published elsewhere [2H]- 

The remainder of this paper is organized as follows. In Section II we describe our 
model system and give the deflnition of DNA-DNA interaction forces. We calculate 
the interaction forces for different counterion and salt concentrations in Section III. 
Section IV is devoted to the free energy perturbation theory for deflning the liquid- 
liquid coexistence densities. We conclude in Section V. 

2. The model 

We consider a B-DNA molecule which is a double helix with a pitch length P=34A. 
There are iVp=20 phosphate charges per one helical turn which makes one elementary 
charge per each 1.7A of axial rise. The geometrical shape and charge distribution of 
DNA is modeled through the accurate Montoro-Abascal model (MAM) [^ Its 
cross section is illustrated in Fig. [H 

A pair of DNA molecules are placed along the xy diagonal of the cubic simulation 
box of volume V = . The size of the simulation box L = 102 A corresponds to 



DNA condensation and redis solution: Interaction between overcharged DNA molecules. 4 



three full turns of DNA. The box also consists of Nq multivalent ions with valency gg, 
= Ns_ + QqNq monovalent colons and = Ng^ + Nc monovalent counterions. 
Here Ns+ = Ng- = Ns is the number of added salt ion pairs, Nc is fixed by the DNA 
phosphate charge due to the constraint of global charge neutrality {Nc = 120 in our 
simulations). All ions are modelled as charged hard spheres. Tetravalent counterions are 
assumed to represent spermine ions. Though the latter is a polyion in vivo, a spherical 
approximation can be used, since the bibliographical data support the idea that it is the 
charge of counterion, rather than its structural specificities, which is important in DNA 
condensation and redissolution processes p. Beside of this, the fluidity of the ordered 
DNA phase also suggests that spermine binds like an ordinary ion to the DNA surface, 
rather than forming inter-strand cross-links to neighboring DNA molecules IT^ . 
Assuming that the ions of the same valency are indistinguishable, the actual number of 
mobile ion species in the simulation box is reduced from five (which are multivalent ions 
and their colons, positive and negative monovalent salt ions and monovalent counterions 
that neutralize the DNA phosphate charge) to three: multivalent counterions (charge 
gg = 4 and diameter dq) and positive and negative small ions (charge Qc = ±1 and 
diameter dc). Periodic boundary conditions in all three directions are applied to reduce 
finite size effects. The DNA replicas in the z direction produce infinitely long molecules. 
The whole system is held at room temperature T = 298K and the water is modelled as 
a continuous dielectric medium with dielectric constant e = 80. A typical snapshot of 
the simulation is shown in Figure |21 

The interaction potentials between the five sort of particles (three of them are the 
mobile ions in solution mentioned above, and two of them, one charged and the other 
neutral, belong to the DNA molecule, see Figure ^ are a combination of hard core and 
Coulomb potentials: 

r oo for r < {di + dj)/2 

forr> (rf, + f/,)/2. 

Here r is the inter-particle distance, i,j= Q for multivalent counterions, c for small ions, 
p for phosphate charges and n for neutral spheres in the MAM (with g„=0). There is 
also an excluded volume potential V^^ between the DNA hard cylinder and the free ions 

hi = Q,c. 

The basic quantity of interest is the effective force per helical turn j^Hl 

F = Fi + #2 + F3. (2) 

between two DNA molecules. Here Fi is the direct Coulomb repulsion per helical turn 
of one DNA molecule as exerted from the phosphate groups of the other DNA. It does 
not depend on salt density and its calculation is straightforward j2Zj- The second term 
F2 in Eq. (j2I) corresponds to Coulomb interactions between the phosphate charges at 
positions ffF {k = 1, ..Np) and the mobile ions at positions r^* (Z = 1, ..Ni, i = c, Q), 

^2 = -ttE ( E E^r^VAl - r7 I)) . (3) 

'^fc=l \ i=c,Q 1=1 I 
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Here < ... > denotes canonical average over small ion configurations. The third term 
F3 in Eq. (0) arises from the entropic contribution of small ions due to their moment 
transform to the DNA surface iS, 

Fs = -hsTj^df (|e^(p.(0))- (4) 

Here / is a surface normal vector pointing outwards the DNA core and pj {j = c, Q) 
is the inhomogeneous microion concentration. The canonical averages in F2 and F3 are 
carried over all configurations of the mobile ions. 

We have performed extensive grand-canonical molecular dynamics (GCMD) 
simulations, similar to those of Ref. ^7\, for a range of different tetravalent counterion 
and monovalent salt concentrations. Each simulation is referred by its nominal 
tetravalent counterion concentration C (salt ion concentration Cg) defined as a ratio 
between the total number of ions Nq (Ns) in the cell without the DNA molecules and 
the system volume V, C = Nq/V (c^ = Ng/V). Additional simulations have been 
carried for these bulk phases in order to calculate the chemical potential fi of solution. 
Then in the simulations with DNA molecules the number of ions in the simulation 
cell was automatically adjusted to the specified value of chemical potential /x. The ion 
diameter was chosen to be dQ=8A for tetravalent counterions and dc=4:A for other free 
ions in the system. This parameter defines the closest approach of the ion to the DNA 
surface and has a strong impact on the polyion electrostatics. The time step At of the 
simulation corresponds to an average ion displacement of O.OSA per time step such that 
the reflection of counterions following the collision with the combined surface of DNA 
is calculated with high precision. About 5 x 10"^ MD steps are required on average to 
reach equilibrium. The time evolution is then followed for 5 x 10^ — 5 x 10^ steps to 
gather statistics to calculate canonical averages of the interaction forces. 

3. Results for the interaction forces 

The case of a single DNA molecule in the presence of spermine ions has been considered 
in Ref. [2H1- It was shown there (see Figure 2 in Ref. that the ionic cloud may not 
only compensate the polyion charge but even exceed it, resulting in an opposite sign of 
the electrostatic potential at some distances from the DNA surface. The adsorption of 
Bjerrum pairs jSHj onto the DNA surface at high tetravalent counterion concentration 
creates consecutive layers of charges of different sign around the DNA molecule. The 
onset of a multilayer structure occurs at C = 1.8mM. Addition of monovalent salt shifts 
this threshold concentration to lower values of C, in accordance with experimental 
observations and two-component Manning condensation theory [211 . For multivalent 
counterion concentrations exceeding C = 1.8mM, the total charge in the DNA grooves 
remains constant and only the total ionic charge adsorbed on the strands contributes 
to the overcharging, similar to our earlier findings [23 I2H]- Beside of this, there is a 
competition between the multivalent and monovalent counterions in binding to the DNA 
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surface as C increases. The multivalent ions tend to replace the monovalent counterions 
on the DNA surface. Thus, at higher C the main contribution to the formation of 
charged layers around DNA molecules comes from Bjerrum association between big 
counterions and small colons. Such charged layers give rise to an attraction between 
two parallel DNA molecules, as shown in Figure El A decrease of the big counterion 
charge leads to the break-up of the Bjerrum counterion-coion pairs and thus destroys 
the layer formation around the DNA molecule. This ultimately results in the loss of 
intermolecular attraction, see simulation data for qq =1,2 and 3 in Figure El 

It should be mentioned that in addition to the intermolecular (or axis-to-axis) 
distance i?, there are angular variables which define the mutual configuration of two 
parallel DNA molecules corresponding to the orientation of their grooves and strands 
[26] . At short surface-to-surface distances between the two DNA molecules up to SA 
there are strong contributions from particular DNA-DNA configurations jSZj. For larger 
separation distance, R > 25A we find no detectable dependence of the interaction forces 
on the mutual orientations of DNA molecules. On the other hand, there is experimental 
evidence that in DNA condensates two neighboring molecules never approach each other 
closer than sA. This apparently means that at such small distances a strong repulsion 
between DNA solvation shells exists. Thus in all figures hereafter we show orientationally 
averaged interaction forces starting from the distance R = 24A. The interaction forces 
are scaled per DNA pitch, i.e. per 10 DNA base pairs. 

The electrostatic F2 and entropic F3 components of total interaction force F for 
tetravalent counterions corresponding to DNA overcharging are separately plotted in 
FigureO It can be seen that the electrostatic force shows oscillations around zero, which 
are reminiscent of the multilayered structure of charges around a single DNA molecule 
[2H1- The deep attractive minimum of the total force F has an entropic origin, whereas 
the second minimum at intermolecular distance R=4lA has a purely electrostatic origin. 

In Figure El we plot the DNA-DNA interaction force for both the undercharged 
and overcharged cases at different salt concentrations. It is evident that in dense salt 
solutions the attractive minimum becomes weaker. For overcharged DNA, C = 65mM, 
the position of the second minimum shifts towards the DNA surface. Thus, whereas the 
DNA overcharging does not strongly depend on the added salt concentration Cg [2Hj, the 
effective forces do. 

The positions of the minimum and the maximum of the force shift towards the DNA 
surface also for higher spermine concentrations. This trend is shown in Figure El where 
the DNA-DNA interaction force is plotted for different spermine and salt concentrations 
at fixed distance R = 38A. For low C, which in Figure El corresponds to the area to the 
left of point A, the DNA-DNA interaction is overall repulsive. Between the points A 
and B a first minimum develops in F. As the spermine concentration increases further, 
the minimum shifts toward the DNA surface and the force has a repulsive tail. This 
tail indicates the appearance of second maximum, which obviously is followed by second 
minimum. 

The dependence of the DNA-DNA interaction force on C for two fixed 
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intermolecular distances are shown in Figure [7| Five different arrows in this picture 
point to different values of C which characterize the number of attractive minima in 
the interaction force F. For small C = O.OlmM, indicated by the arrow next to a 
in Figure [71 the interaction force has no minimum and thus is totally repulsive. For 
spermine concentration C = 0.025mM, corresponding to the arrow next to b, the force 
has a single minimum. For intermediate C = 1.7mM and C = 56mM, see the arrows 
next to c and d respectively, there are two minima in the interaction force (note that a 
positive force at i? = 38A for C corresponding to arrow d indicates the occurrence of 
second maximum, which obviously if followed by second minimum at larger R). And 
finally, at even higher C = 400mM, the first minimum has disappeared, however the 
second minimum is retained. 

A full set of distance resolved DNA-DNA interaction force curves for different 
tetravalent counterion densities C are presented in Figure |H1 It can be seen that even 
a small trace of spermine ions - well below the overcharging threshold - induces an 
attraction between the DNA molecules, except at very close distances, see curve for 
C=0.1mM. As the DNA molecules get more overcharged, which corresponds to high 
spermine concentrations, the main minimum narrows and becomes more shallow. At 
the same time, the width and the height of the maximum increases. We note that the 
attractive minimum for undercharged DNA pair has mainly a pure electrostatic origin 
and arises due to charge correlations in the electrolyte. However, for overcharged DNA 
the main contribution to the force at this minimum is due to the spermine crowding 
near the DNA surface. The second maximum originates both from electrostatic and 
entropic forces. Finally, the second minimum emerges from a pure electrostatic effect. 

In the following we calculate the total effective pair interaction potential U{R) 
per unit length for a given bulk salt concentration Cg and different Spe concentrations 
C. The quantity U{R) is obtained by integrating the distance-resolved interaction 
force averaged over all microion configurations [HI]. Results are shown in Figure IHl 
The oscillations in the force imply that the interaction potential also oscillates. With 
increasing C, the first minimum of U{R) is getting deeper and is achieving a maximal 
depth at the overcharging concentration C ~ 1.8mM. A further increase of C again 
reduces the depth of this minimum. The position of the minimum, on the other hand, 
hardly depends on C. The Spe-layering around the pair of DNA molecules induces a 
second minimum at larger separations as revealed in the enlarging inset of Figure El This 
minimum is of electrostatic origin and occurs for C >65mM. Again the depth of the 
second minimum first increases and then decreases with increasing C. At intermediate 
Spe concentrations, we are thus confronted with a double minimum potential which is 
induced by layering. It is worth to mention that positions of the second minimum can 
be related to intermolecular distances between the DNA molecules in cholesteric phase 
induced by polyamines j32l 0] • Direct measurements and theoretical investigations 
PT] of intermolecular forces demonstrated that the energetics of this cholesteric phase 
is determined primarily by electrostatic interactions. Since the potentials in Figure IHl 
are scaled for one DNA pitch length, the interaction strength corresponding to the 
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minimum of curve (2) for very long DNA molecules is sufficiently large to induce 
condensation. This implies that DNA aggregation can take place well below overcharging 
Spe concentrations. 




dR. (5) 



4. Phase diagram for double-minimum potential. 

The characteristic double-minimum structure of the interaction potential U{R) will give 
rise to an unusual phase behavior. We have calculated the phase diagram of a columnar 
DNA assembly on the basis of our simulated effective pair interactions. To do so, we 
assume that the DNA molecules are parallel along a certain length ^. This length is an 
additional parameter which we fix to be £ = 20 x P. We comment on the dependence of 
the phase diagram on ^ in Section V. The assembly of parallel DNA can be considered 
as a two-dimensional many-body system interacting via a potential £ x U{R) and being 
characterized by a DNA particle number density p. We calculated the free energies of 
the fiuid and solid phases by using different techniques outlined below and perform the 
traditional Maxwell double tangent construction to identify the coexisting densities. 

The free energy of dilute fluid phase is approximated by the two-dimensional 
perturbation theory [33] via splitting the total potential into a repulsive and attractive 
parts, U{R) = Ur{R) + Ua{R)- The repulsive potential Ur{R), identical to U{R) but 
truncated and shifted towards zero at the first minimum at i? = Rmin, is mapped onto 
that of effective hard disks of diameter 

cRri , 

1 — PTfn I — 

Here the cross-section diameter for the DNA molecule is a = 20A. The total Helmholtz 

2 

free energy involves that of a hard disk fiuid with effective area fraction rj = — and 
a mean-field correction which we simply model as vrp 

joc u^^^^_ Yoi the free energy 
of a hard-disc fiuid, analytical expressions are available [37J. The free energy of the 
solid phase, on the other hand, is calculated as a lattice sum with the assumption of a 
two-dimensional triangular lattice structure. The lattice constant is directly related to 
the DNA number density p. 

Figure ITUl shows the resulting phase diagram with the coexisting DNA densities for 
the wide range of C . At low C there is a strong first-order gas-crystal phase transition 
whichs coexistence region is widened for larger C due to the increasing attractions. 
Between the two threshold concentrations C ~ O.lmM and C ~ 65mM there is enough 
attraction to stabilize a liquid phase of high DNA density. At higher Spe concentrations 
a second crystal, with a considerable larger lattice constant than that of the high-density 
solid emerges. We call this novel phase a mesocrystal since its density is intermediate 
between that of the fiuid and the almost closed-packed solid. 

Condensation (see the cross in Figure ITUj) and subsequent redissolution (dot-dashed 
line in Figure EH) of dilute DNA solution are other implications of phase diagram. As the 
spermine concentration increases for fixed pa^ = 0.002, which corresponds to a typical 
DNA concentration of Img/ml DNA, first the gas-liquid coexistence line is encountered. 
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This is the condensation transition and the system will split into a low density gas and 
a high density liquid phase. At much higher C the coexistence line is touched again 
and the system redissolves back into the dilute gas phase. The corresponding threshold 
concentrations of the condensation and redissolution are in the range Cc ~ 0.3mM and 
Cd ~ 165mM and agree well with the experimental observations [?H ITT]. 

5. Discussions and conclusions 

One issue we want to discuss is the dependence of the phase diagram on the DNA length 
^. Since £ is a prefactor of the effective potential, it plays formally the role of an inverse 
system temperature. We have explored the phase behavior for smaller and larger DNA 
segment lengths of £ = 5P and £ = lOOP respectively. Firstly, the stability of liquid 
pocket depends sensitively on £; it disappears completely for small but extends towards 
larger C for larger ^. Secondly, the fluid coexistence density shifts to considerably higher 
values for smaller C.. Hence, condensation and redissolution is prohibited for small DNA- 
segment lengths. This is in line with the experiments of Ref . jHH] , where a threshold value 
of £ ~ 15P for the minimal length ^ required for condensation is reported. For a triplex 
DNA (three stranded DNA molecule) the minimal length is reported to be about 2P 
jSH] • The difference between the minimal lengths for duplex and triplex DNA segments 
arises from the DNA charge density [311 EDI- The higher the linear charge density, the 
stronger the overcharging. This in turn results in the precipitation of triplex DNA at 
spermine concentrations, where duplex DNA does not aggregate |3Hl 112] • 

The DNA fragments considered here have a length less than the persistence length 
of the DNA molecule which is around 500A. That is why one can safely adopt a rigid 
rod assumption for the DNA molecule and avoid the intramolecular fluctuations (off 
the long DNA axis). This is also in line with the experiments [11-16] where the DNA 
redissolution measurements were done for DNA segments which are smaller than the 
persistence length. Hence the particles involved are very stiff rods. 

We speculate that the height of the energetical barrier between the two minima in 
the interaction potential U{R) could regulate the DNA segment lengths in the crystalline 
structures observed in experiments. Imagine that the solution consists of a mixture of 
DNA segments of different lengths (but still smaller than the persistence length 500A). 
At lower spermine concentration, when the interaction potential has a single minimum, 
all DNA segments will form a bundle with an average inter-particle distance of the order 
of R = 28A. However, when the energetical barrier develops at intermediate spermine 
concentrations with a height around O.bksT — Ifc^T per pinch length, only short DNA 
segment could overcome this repulsion. Thus, short segments will fall into first minimum 
and form a dense hexagonal structure, while longer DNA segments will be trapped in the 
second minimum and form more swollen, fluid-like structure, apparently a cholesteric 
phase. 

The other issue is the orientational entropy of DNA molecules and its influence 
on the interaction potential and free energy of DNA solutions. The evaluation of the 
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contribution to the free energy from the oricntational entropy of two nonparallel DNA 
rods (intermolecular fluctuations) in MD simulations is a tremendous task. In fact, the 
DNA-DNA effective potential must be averaged over all possible orientations for two 
DNA. Unfortunately, an implementation of tilted DNA molecules will break the system 
symmetry in simulations. On the other hand, simulation of 20P long molecular segments 
for the system parameters invoked here still is a challenge. For the bundle phase, due to 
the deep attraction (more than 40 ksT) between the DNA rods of 20P length, all other 
than parallel configurations of DNA molecules have little statistical weight. Thus, only 
a tiny correction to the energy of bundle is expected. 

For the dilute gas phase the oricntational entropy for DNA rod has an upper limit 
of about 2k bT. At lower DNA densities, when the average DNA-DNA distance is of the 
order or less than the DNA length, this entropy should be added to the free energy of 
liquid phase. On the length-scale of free energy, where both the liquid perturbation and 
lattice sum produced curves have deep minima about dozen fc^T at intermediate DNA 
densities, such a upward shift of hquid free energy at smaller DNA densities will have 
only a slight effect on the phase diagram for 20P long DNA fragments. In other words, 
a corresponding upward shift of the free energy of dilute DNA system will not strongly 
affect the coexistence spermine concentration values deduced from the double-tangent 
procedure. Hence, a length of 20P for the DNA segment is enough to safely neglect 
oricntational entropy effects in the phase diagram shown in Fig. 10. 

In conclusion, we have calculated the influence of tetravalent countcrions on 
the effective interactions and the phase diagram of columnar DNA assemblies by 
primitive-model-type computer simulations and statistical theories. We find that 
a small concentration of tetravalent counterions induces DNA condensation. The 
layering of the strongly coupled tetravalent counterions on the DNA strands yields 
an oscillatory effective interaction potential with a double-minimum structure at 
intermediate counterion concentrations. This explains the redissolution transition and 
triggers a novel stable mesosolid. Our threshold concentrations are in good agreement 
with experimental data. 
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Figure captions 



Figure 1. (Color online) Cross section of DNA in the xy plane for the Montoro- 
Abascal model (MAM). Phosphate charges are shown as dark spheres. The DNA 
cylindrical core is colored in gray, the hatched areas correspond to neutral hard 
spheres. The inscribed letters "M" and "m" denote the major and minor grooves 
correspondingly. 



Figure 2. (Color online) Typical snapshot in the simulation cell. The DNA molecules 
are shown as two parallel rods in the z direction, over-wrapped by two strings of light 
grey (neutral sphere in MAM (see text), colored yellow in online figure) and grey 
(phosphate sphere in MAM, red in online figure) spheres. The tetravalent Spe ions are 
shown as big black (blue in online figure) spheres. Light grey (yellow in online figure) 
spheres represent colons, and dark grey (green in online figure) spheres are monovalent 
counterions. 



Figure 3. (Color online) DNA-DNA interaction force F/Fq versus intermolecular 
separation distance R for Cs=25 mM and C=65 mM. The charge of the big counterions 
is indicated next to the corresponding curves. Fq = ksT/P, where P = 34A. 



Figure 4. (Color online) DNA-DNA interaction force F/Fq versus intermolecular 
separation distance R for Cs=25 mM, C=160 mM. Fq = ksT/P, where P = 34A. The 
parameters corresponds to complete DNA overcharging. Note that there is a second 
maximum in the total force F at about R=36 A and a second minimum at about 
E=4lA. 



Figure 5. (Color online) DNA-DNA interaction force versus intermolecular separation 
distance R for ^,=25 mM (black lines in online figure) and Cs~320 mM (red lines in 
online figure) and two different spermine concentrations: C=0.8mM (imdercharged 
DNA, solid hues) and C=65mM (overcharged DNA, dashed hues). Fq = ksT/P, 
where P = 34A. The curves that correspond to C=65mM are shifter upward for 
clarity. 



Figure 6. (Color online) DNA-DNA interaction force versus spermine concentration 
C. The intermolecular separation distance is i?=38A (roughly the position of the 
second maximum for C=65mM). In the region between the points A and B (for salt 
concentration Cs = 100mM) first minimum of the interaction force develops. Fq = 
ksT/P, where P = 34A. 
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Figure 7. (Color online) DNA-DNA interaction force versus spermine concentration 
C at fixed intermolecular separation distances R for Cs=25mM. At R =30A (R = 38A) 
the first (second) minimum emerges at intermediate spermine concentrations C and 
then disappears at higher C. Arrows and labels from a to e next to them are a guide to 
eyes to point the different spermine concentrations where the shape of the interaction 
force undergoes significant changes. For more details see the text. 



Figure 8. (Color online) DNA-DNA interaction force versus intermolecular 
separation distance R for Cs=25mM, Different spermine concentrations are indicated 
next to corresponding curves: C=0.1mM, 0.8mM, 18mM, 65mM, 160mM, 280mM, 
400mM. Fo = ksT/P, where P = 34A is the DNA pitch length. 



Figure 9. (Color online) DNA-DNA effective pair potential versus intermolecular 
separation distance R for Cs=25mM. The spermine concentrations arc C=OmM (1), 
O.lmM (2), 0.8mM (3), 18mM (4), 65mM (5), 160mM (6), 280mM (7), 400mM (8). 



Figure 10. (Color online) Coexisting DNA densities as a function of Spe 
concentrations C, for Cs=25mM. The stable phases found are gas-like, liquid and two 
triangular crystals with different lattices constants (solid I and solid II). All phase 
transition between these phases are of first order. For the sake of better resolution at 
smaller DNA densities, the 2^axis is expanded below pa"^ < 0.015. 



